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Abstract 

Quantum field theory of a damped vibrating string as the simplest dis- 
sipative scalar field investigated by its coupling with an infinit number 
of Klein-Gordon fields as the environment by introducing a minimal 
coupling method. Heisenberg equation containing a dissipative term 
proportional to velocity obtained for a special choice of coupling func- 
tion and quantum dynamics for such a dissipative system investigated. 
Some kinematical relations calculated by tracing out the environment 
degrees of freedom. The rate of energy flowing between the system 
and it's environment obtained. 



1 Introduction 

In classical mechanics dissipation can be taken into account by introducing a 
velocity dependent damping term into the equation of motion. Such an ap- 
proach is no longer possible in quantum mechanics where a time-independent 
Hamiltonian implies energy conservation and accordingly we can not find a 
unitary time evolution operator for both states and observable quantities 
consistently. 
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To investigate the quantum mechanical description of dissipating systems, 
there are some treatments, one can consider the interaction between two sys- 
tems via an irreversible energy flow [1,2], or take a phenomenological treat- 
ment for a time dependent Hamiltonian which describes damped oscillations, 
here we can refer the interested reader to Caldirola-Kanai Hamiltonian for a 
damped harmonic oscillator [3]. 

H [t) = e~ m ^- + e 2 ^-mu 2 q 2 (1) 
w 2m 2 y w 

There are significant difficulties about the quantum mechanical solutions of 
the Caldirola-Kanai Hamiltonian, for example quantizing with this Hamil- 
tonian violates the uncertainty relations or canonical commutation rules [4] 
and the uncertainty relations vanish as time tends to infinity and this is be- 
cause the related time evolution operator is not unitary. 
In 1931, Bateman [5] presented the mirror- image Hamiltonian which consists 
of miror Hamiltonians, one of them represents the main one-dimensional 
damped harmonic oscillator. Energy dissipated by the main oscillator com- 
pletely will be absorbed by the other oscillator and thus the energy of the 
total system is conserved. Bateman hamiltonian is given by 

H = — + -^—(xp — xp) + (k — — )xx, (2) 
m 2m Am 

with the corresponding Lagrangian 

L = mxx + ^-(xx — xx) — kxx, (3) 

canonical momenta for this dual system can be obtained from this Lagrangian 

as 

dL ± 0_ dL . (3 

p = — = mx - -x, p = — = mx + -x, (4) 

ox 2 ox 2 

dynamical variables x, p and p, x shoud satisfy the commutation relations 

[x,p] = i, [x,p]=i, (5) 

however the time-dependent uncertainty products obtained in this way, van- 
ishes as time tends to infinity. 
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Caldirola [3,6] developed a generalized quantum theory of a linear dissipa- 
tive system in 1941 : equation of motion of a single particle subjected to a 
generalized non conservative force Q can be written as 

d f dT^ dT dV , 

where Q r = —/3(t) J2 a rjQj, and a rj 's are some constants, changing the vari- 
able t to t* using the following nonlinear transformation 

t* = x (t), dt = 4>(t)df, ( p(t) = Jl^ dt \ (7) 

and from the definitions 

rJn r)T* 
? = ^ L* = L( q ,f,n P*=Qf> (8) 

the Lagragian equations can be obtained from 

op) ~ -q^ = °' ( 9 ) 

where H* = J2p*Q* — L*. Canonical commutation rule and Schrodinger 
equation in this formalism are 

[q,p*]=i, H*iP = t^, (10) 

but unfortunately uncertainty relations vanish as time goes to infinity. 
Perhaps one of the effective approaches in quantum mechanics of dissipa- 
tive systems is the idea of considering an environment coupled to the main 
system and doing calculations for the total system but at last for obtaining 
observables related to the main system, the environment degrees of freedom 
must be eliminated. The interested reader is referred to the Caldeira-Legget 
model [7,8]. In this model the dissipative system is coupled with an envi- 
ronment made by a collection of iV harmonic oscillators with masses m n and 
frequencies u> n , the interaction term in Hamiltonian is as follows 

N N 2 

where q and x n denote coordinates of system and environment respectively 
and the constants c n are called coupling constants. 
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The above coupling is not suitable for dissipative systems containing a dis- 
sipation term proportional to velocity. In fact with above coupling we can 
not obtain Heisenberg equation like q + u 2 q + (3q = £(£) for a damped har- 
monic oscillator consistently and we can not study dissipative quantum fields, 
for example, a dissipative vibrating medium with this model. In this paper 
we generalize the Caldeira-Legget model to an environment with continuous 
degrees of freedom by a coupling similar to the coupling between a charged 
particle and the electromagnetic field known as minimal coupling. In sections 
2, the idea of a minimal coupling is introduced and the quantum dynamics 
of a damped vibrating string as the simplest scalar field theory, is investi- 
gated. In section 3, quantum dynamics of the string and it's environment 
is investigated. In section 4 some transition probabilities indicating the way 
dissipation flows, are obtained. 



2 Quantum dynamics of a damped vibrating 
string 

In this section we consider a damped vibrating string as the dissipative sys- 
tem although the method is general and can be applied to a general scalar 
field. Quantum mechanics of a damped vibrating string with mass density 
A, tention jj, and lengh L, can be investigated by introducing a reservoir or 
an environment that interacts with the string through a new kind of mini- 
mal coupling. Let the two ends of the string be fixed in x = and x = L 
respectively and vibration be only in the y direction. If ip(x, t) is the wave 
function of the string, to quantizing ift(x, t), we assume ift(x, i) to be a hermi- 
tian operator and can be expanded in terms of orthogonal functions, sin 
as follows 

Mx, t) = E -^=[a n {t) + a{{t)\ sin (12) 

where u n = \J\^ and a n and a} a are annihilation and creation operator 
of the string respectively and satisfy in any instante of time the following 
commutation rule 

[a n (t),at n (t)] = 6 nm , (13) 
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By definition of a conjugate canonical momentum density as 



n f (x, t)=ij>2\J - a n (t)} sin (14) 

then from ()13|) ip, tc^ satisfy equal time comutation relation 

[if)(x, t), tt^(x', t)] = z8(x — x'). (15) 

Hamiltonian of string is defined by 

r L tt 2 i 1 00 1 

H s = J o dx(— + -fMipl) = E ^K-a n + ^) ( 16 ) 

Let the total Hamiultonian, i.e., string plus environment be like this 

rT / n f L , (Tc^ix.t) - R(x,t)) 2 1 , 9 

= ^ dx y ^ 1 + + Hb) (17) 

where ^ denots derivative with respect to x and ^ is a constant depending 
on string properties, Hb is the reservoir Hamiltonian 



00 r+°° 1 1 

= E / fto&l&Wdt) + 5), w s = |*|. 



;is) 



Annihilation and creation operators 6 n g, & in any instant of time, satisfy 
the following commutation relations 

[Kk{^l%M = 5n m 8{k-k'), (19) 

and we will show later in section 3 that reservoir is an infinit number of 
independent Klein-Gordon equations with a source term. Operator R(x, t) 
have the basic role in interaction between string and reservoir and is defined 
by 

r+00 , T17VT 

R(x, f) = £ / d 3 k[fMKk(t) + fMtfjt)} sin —, (20) 

n=1 J-oc Li 

let us call the function f(wk), the coupling function. By using of (|TH|) . it can 
be shown easily that Heisenberg equation for ip(x, t) and 7r$(x, t) leads to 

^(x,t)=i[H,^(x,t)] = 7 ^j-^, 

7r^(x, t) = i[H, 7iip(x, t)) = nii xx , (21) 
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which after omitting 7ty, gives the following equation for the damped vibrat- 
ing string 

Xi> - fM^xx = -R(x,t). (22) 
Using (JT5)l the Heisenberg equation for b n ^ is 

, . _ ;[ff h J — - _j_ ;f*(,.,J\ I ,/,/V A ^-n _ 

nfc 



777TT 

• 6 J = -iwjft^ + jf sin —da/, (23) 



with the following formal solution 

&njf(*) = 6 njf(°) e " <WSt + i /*K) tdt'e-^-n [ L ^ x ',t') S in '-^dx', (24) 



mix' 



substituting 6 n ^(t) from (J2IJ) into ((22), using the relation 8{x—x!) = \ E£°=i sin sin ^ 
and at last integrating over x' gives 

A^ - ^ + f dt'ip(x, t)i(t - f) = £(x, t) , 
Jo 

^(x, t) = il <Pku>&(uz)bJp)e-*# - fK)^(0)e^] sin — 



oo 
oo 



(•OO 

7(t)=47rL/ dwjfl/^l^lcoswjft, (25) 

w 

it is clear that the expectation value of £(x,t) in any eigenstate of Hb, is 
zero. For the following special choice of coupling function 



\ 



^ (26) 



4n 2 Lu, 

k 



equation (|23j) takes the form 

A^ - fi^ xx + f3ip = i(x, t), 

«(«. o = i/5r £T S |6 « j(o)e "" s ' - "-s* '^' 1 81,1 (27) 



Hisenberg equation for a n and a* is 



V ' L\uj n V LA n A 7- 

• Xu) n ..± . 2 (a) + a n ) 
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Definition A n = ^0= and B n = i\f^(al - a n ) and 



easily obtain 



using 



we can 



A n + u 2 n A n + ^A n = Cn(t), 



Cn(t) = i\ 



(3 r+™ d 3 k 



4vr 2 A 2 L 
with the following solution 



"s* - &^(0)e^] 



(29) 



Mt) = e-£{E n e iQnt + F n e- lfW ) + M{t), 

"r U°) 



M n (t) = i / drk 



+ OO 



/3 



\47r 2 A 2 L^ 



(30) 



where f2 n = ycu 2 — Operators i£ n and F n , are specified by initial condi- 
tions 



E n + F n = A n (0)-M n (0), 
(Z^ + ^ n )E n -(^+^n)^ = A»(0)-M n (0), (31) 

solving above equations and substituting E n and F n in (|5UJl one obtains 
A n (t) = e-^{A n {0) cosQj + Jl-^O) sinfi n t - ^p s inQ n i 

-M n (0) COSn n t + A n(U)-M n (0) sinnj y + ^2) 

i in 



It is clear from (|28jl that A n (0) is dependent on string and reservoir operators 
in t — . Substituting A n (t) from (}3*2"|) in (|2*Hl we can find a stable solution 
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for b „j?(t) in t — > oo as 



+iu^(A n (0) - M„(0))} 



L/3 



We 



■KA(0)-^M n (0) 



i/3 



V A: 



+00 



6^(0) sin^^t l^rf 

e } • 



(33) 



now substituting 6 from in (j^Sjl and using (J3~2j) . one can show that 

OO r+oo 

B n (t) = \A n (t) + £ / d 3 k^ 



n=l ' 



\ 



/3 



Ati 2 LujI ' 



(34) 



A vector in fock space of string can be written like this 



00 00 



l $ )* = E E $ni,..,«>l), 

j=0 m,...,n,j=l 



\ n j) 



lis 



(35) 



where \n) a denotes the state of a single particle in mode n, with corresponding 
wave function (x\n) = y^sin 1 ^. Operators a n (0) and 0^(0) act on basis 
vectors \ni) s <g> ... ® |7i 3 -) a = |ni, ...,rij) s as follows 



a m (0)|ni, ...,nj) a = 5nr,m\ni, ...,n r _i,n r+ i, ...,n,) s , 

r=l 

a^(0)|m, ...,n,) s = |m,m, ...,%)„, (36) 
also a vector in fock space of reservoir can be written as 



00 00 



^>B = E E / ^•••^*n 1) ... ) n j A,--- ! fc J )|^l ) ^l)B® 



j=Q vi,...,i/j=l 



(37) 

In subsequent section we show that reservoir is infinit number of indepen- 
dent Klein-Gordon fields and we can interpret \k,v)s as a single particle 
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state belong to uth Klein-Gordon field with corresponding momentum, k. 
Operators b n ^(0) and 6^(0) act on basis vectors \ki,Vi)B ® ••• <£> \kj,Vj)B = 
|A?i,z/i, i/j-Jb as 



6 n f(0)|A;i, z^i, z^)b 

1 ! 7 ■• • j kj , Uj ) b ) 

r=l 

^0)1^1) "i> = W,n',k 1 ,u ll ...,k jl u j )B- (38) 



If the state of system in t = is taken to be | -0(0) ) = |0)b <8> ■■■m r ) s 
where |0)b is vacuum state of reservoir and \rrii, ...m r ) s an excited state of 
the Hamiltonian H s then by making use of ([32)1 . ()33|) and (|34|) it can be 
shown that 

f L 1 ■ 1 

fim^-Kx,[B(0| ® s (mi, ...,m r | : / dx[-\ip 2 + ~/^] : l m i> •■•> m r) s ® |0)b] = 0, 



o 



^m^oo[B(0| (g) s (mi, ...,m r | : ^ u n al(t)a n (t) : |mi, ...,m r ) s ® \0) B ] 



n=l 

L[3 2 uJn . f+°° dx e ixt 



r 

x,(m r , ...mi| : A^(0) : |m 1; ...m r ) s } ~ £ — . (39) 



2 167T 2 A i-oo X Co" 2 — X 2 + 

/3 2 1 



8A i=1 o; mi 



where : : denotes the normal ordering operator. Now by substituting b At) 
from into (JTHJ), we have 

limt-^oo[B(0\ <S> s (m r , m\\\ J d 3 k^ u k b ik b nk '■ l m i> ® |0) B ] 

/3 3 roc dx 



71=1 



1=1 • yu 



2ttA mi -A) - x 2 ) 2 + gx 2 

4 \ . / 



r /"CO r v^A r r 

, 2 . (40) 

2vrA^ mi io « - x 2 ) 2 + g* 2 
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3 Quantum field of reservoir 

Let us define the operators Y n (x,t) and H n (x, t) as follows 



/+oo d^h — 
, {b nn {ty k£ + b\(t)e~ ik ^ 

+00 



!I„(.F. /) = / / d"kJ 1 ^{b[ } .{t)e-^ - bjfte**), (41) 



CX) 



then using commutation relations ()19j) . one can show that Y^(x, t) and Il m (x, £), 
satisfy the equal time commutation relations 

[Y n (x, t),Tl m (x>, t)] = i5 nm 5(x - x'), (42) 

furthermore by substituting b %(t) from ( l2*4*j) in (}4~T|) we obtain 

= v 2 r n + Li n (t)P(f), P(f) = Re £J d'kJ^^fi^e-^, 

BY ■ r+°° f(uA 

U n (x,t) = - LA n (t)Q(x), Q(x) = Jm / rf 3 fc M fcJ e"^, 
^ ^/2(2vr) 3 ^ 

(43) 

so for any n, l^(x, t) satisfies the following source included Klein-Gordon 
equation 

^ - V 2 F n = LA n (t)Q(x) + 2A n (t)P(x), (44) 
ot z 

with the corresponding Lagrangian density as follows 

£n = 2 ( ^f )2 " 2 VF "- Vy " " LA nQW^f + LA n P(x)Y n . (45) 
It is clear that the reservoir is made by an infinit number of massless Klein- 
Gordon fields containing the source term 2A n Q(x) + 2A n P(x). Hamiltonian 
density for (|4*4*jl is 

K = (n " + ^ n9)2 + \\VY n \ 2 - LA n P(x,)Y n , (46) 

and equations (|4*3j) are Heisenberg equations for fields Y n and Il n . If we 
obtain b n ^ and 6^- from (}4"T|) in terms of Y n and Il n and substitute them in 
(HBJ), welind 

1 1 TT 2 1 

Hb = J d 3 k^b nR + -) = -f + -IVKI 2 . (47) 
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4 Transition probabilities 

We can write the Hamiltonian (|20|) as 
H = H + H', 

oo y 00 roo I 

H = H S + H B = J2 (4«n + :K + £ / d3k ^k(KAk + 9 ) 

n=l 2 n=l J -°° Z 

H> = -f L dx?^R{x,t) + (48) 
jo A 2A 

and in interaction picture we can write 

Mt) = e iHot a n (0)e- iHot = a n {Q)e~^\ 

Ku(t) = e^b^e-^ = brf(0)e-** t (49) 

terms jir^, and ^ are of the first order and second order of damping re- 
spectively, therefore for sufficiently weak damping, Trr is small in comparison 
with y 7r^. Furthermore |v has not any role in those transition probabilities 
where initial and final states of Hamiltonian of vibrating string are different, 
hence we can neglect the term and estimate if' by — y7i>. Substituting 
a n j and from (}4*9~)) into — j7ty, one obtains if} in interaction picture, as 



+/>s)«l (0)6^(0)e*^«)* - /(o; fc )a„(0)6^(0)e-K+-n)* 
-r(^)a n (0)^(0)e i K~^), (50) 



terms containing just «n(0)6 n ^(0) and ajj(0)6' -(0) violate the conservation of 
energy in the first order perturbation, because a n (0)6 n g(0) destroys an excited 
state of string while at the same time destroying a reservoir excitation state 
and a}j(0)&' -(0) creates an excited state of vibrating string, while creating 
an excited reservoir state at the same time, therefore we neglect the terms 
involving a n (0)6 n ^(0) and ajj(0)6' -(0), because of energy conservation and 
estimate H\ by 

j It 00 r+°° 
-00 



1 / / r+00 
Hj = £ ^^[/K)4(0)^(0)e l 

-r(^)a n (0)^ (Oe"^-^]. (51) 



n=l ' 
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The time evolution of density operator in interaction picture is as follows [10] 

Pl (t) = f/ / (t,t )p / (t )^ t (t,to), (52) 
where Ui is the time evolution operator, which in first order perturbation is 

rt 
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U I (t,t = 0) = l-i f db 1 H' I {t 1 ) = 
Jo 

1 - W A loo A Ev^[/M«i(o)yo)e^- 



n=l 



-fK)q B (Q)^(0)e^] (Mw _; g) . (53) 

2 

Let p/(0) = \m,...,m) r s r s (m,...,m\ <S> |0)b b(0| where \0)b is the vacuum 
state of reservoir and \m, m) r s is an excited state of vibrating string, from 
now on by \m, m) r s , we mean a string state containing r phonons of mode 
m, substituting £//(£, 0) from (|5H|) in (JH2|) and taking trace over reservoir 
parameters, we obtain 

p sI (t) := Tr B (pi(t)) = \m,...,m) r s r s (m,...,m\ 
Lu r+™ sin 2 ^ Wm) t 

(54) 



where we have used the formula Tr B [\k, n) b b(&', s|] = S ns 5(k — k'). In large 

sin 2 

time approximation, we can write — = 2Txt5{u^ — u) m ), which leads 

(-2-2 — ) 2 

to the following relation for density matrix 

U\ I \r r I I , ^^ 7r2u rn^\f( u, m)\ 2 , Nr _i r _i, , 

Paj(t) = |m, ...,m) s s {m, ...,m|H \m, ...,m) s s (m,...,m|, 

(55) 

from density matrix we can calculate the probability of transition \m, m) r s — > 
|m, m) r ~ l as 

>r- 1 r— 1, 



r K...,m>;H'n,...,m>r 1 =Tr [l m '---' m )« « (m, m|p(t)] 

A 

A 



iv.[K... > m>r 1 rv,-,™^)] = 2L ^\jm\\ (56) 
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where Tr s means taking trace over string eigenstates. For the special choice 
(|26|) . above transition probability becomes 

T \m,...,my a ^\m,...,m) r s -i = ^! (57) 

which shows that the rate of phonon number reduction (energy flow), is con- 
stant. Now consider the case where the reservoir is an excited state in t = 
for example p/(0) = \m,...,m) r s r r s (m, m\®\pi, v x , ...pj, v 3 ) B B {pi, v x , —Pj, v. 
where \p\, v\...pj, Uj) B denotes a reservoir state containing j quanta with mo- 
menta pi,...,pj belonging to the z/i,z/ 2 ,...,^th, Klein-Gordon field, respec- 
tively then by making use of 

Tr B [bl£\pL,v 1 ...p j ,u j ) B B{Pi,Vi-Pj,Vj\b m g) = 5 nm S(k-k'), 

j 

Tr B[b n j:\pi,vi-Pj,Vj)B B (j)i,Vi...pj,Vj\$ mj; ] = J2 S n,»i S ™,v l 5(k - pi)5(k' -{ 

1=1 

we find 

Psi(t) = \m, ...,m) r s r s {m, ...,m| + 

L j sin 2 ^prz^rlf 

+ -^J2^rWr,m,...,m) r s r s {u r ,m,...,m\\f(uj fr )\ 2 — 



j i 



2 



r+°° sin — — — / 

+ —\m,...,m) r - 1 r J 1 (m,...,m\ J d*k\f(^)\ 2 —^—^59) 

which gives the transition probability for \m, ...,m) r s \m, ...,m)^ _1 and 
\m, ...,m) r s — > \p, m, m) r s , respectively, as follows 

r K ..., m> ^ K ..., m> - = I>.k-.,m)r 1 :-\m,...,m\ PsI (t)] = ^^\f( 
^\ m ,...,my 3 ^\u, m ,...,my s = Tr a [\v,m, ...,m) r s r s {\v, m, m\p sI {t)\ 



2A p=1 



For the choice we find 



r . 

|m,...,m)J— »|m,...,m) J 1 2A' 

r[ m) ..., m );_»|i. )m ,..., m )r = ^— ^- ^ ^(^Pr ~~ u v)- (61) 

07TA , r=1 UjJ r 
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Another important case is when the reservoir has a Maxwell-Boltzman dis- 
tribution, so let p/(0) = | m, m) r s r s (m, m\ <S> p\ where 

Pb = — ' — -u B , then by making use of following relations 

TR B {e-F?T) 

TtbKwIK^] = TrB^plbl^ = 0, 

Tr B [b nj ;p T B bl g } = S nm 6(k - k')e~^ , 

Tr B {b ] n%P T B b m] :] = S nm 5(k - k'), (62) 

we can obtain the density operator p s i(t) in interaction picture as 

p sI (t) := Tr B [pi{t)\ = |m, ...,m) r s r s (m, ...,m\ 

L 00 [+°° . sin 2 Kr"")^ 

+jvY, UJ n\n,m,...,m) r s r s (n,m, ...,m\ / d 3 k\f(^) | 2 — — £— e «t 
4A n=i ) 

Llo r+°° sin 2 (LJg ~" m) t 

°° I 2 / 

which accordingly gives the following transition probabilities in long time 
approximation 

r \m,..., m )r-+\ m ,..., m y-i =Tr s {\m,...,m) r ~ 1 ^(m, m\p sI {t)\ = ^^ l/k)! 2 , 
r\m,..., m y a ->\u,m,..., m y a = Tr s [\u,m, ...,m) r s r s {v, m, m\p sI {t)\ = 



2Ln 2 uit 



" l |/(^)| 2 e"^, (64) 



A 

substituting (J2fij) in these recent relations, we find 

fit 

2A 

fit 



T ' _ 

|m,...,m>^|m,...,m)s _1 ~~ 



r| m ,..., m >j- > |,y,m,...,m>j = «T. (65) 

So in very low temperatures the energy flows from oscillator to the reservoir 
by the rate and no energy flows from reservoir to oscillator. 
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5 Concluding remarks 



The Caldeira-Legget model generalized to the case where the environment 
has continuous degrees of freedom, for example, a Klein-Gordon field or an 
infinit number of Klein-Gordon fields. A minimal coupling method intro- 
duced which leads to a consistent investigation of the quantum dynamics of 
a large class of quantum dissipative systems. By choosing different coupling 
functions in ([20)1. we could investigate another forms of dissipation. The rate 
of energy dissipation (energy flowing between the system and it's environ- 
ment), was a constant. This problem can be extended to the case where the 
field R, becomes a general field for example a vector field, which is suitable 
for investigating three-dimensional quantum dissipative models. 

References 

[1] H. Haken, Rev. Mod. phys. 47 (1975) 67. 

[2] G. Nicolis, I. Prigogine, Self-Organization in Non-Equilibirium system, 
Wiley, new York, (1977). 

[3] P. Caldirola, Nuovo Cimento 18 (1941) 393. 

[4] I. R. Svinin , Teor. Mat. Fiz. 27 (1972) 2037. 

[5] H. Bateman, phys. Rev. 38 (1931) 815. 

[6] P. Caldirola, Nuovo Cimento 77 (1983) 241. 

[7] A. O. Caldeira, A. J. Legget, Phys. Rev. Lett. 46, (1981)211. 

[8] A. O. Caldeira, A. J. Legget. Ann, phys.( N.Y.) 149,(1983) 374. 

[9] W. P. Schleich, Quantum optics in phase space, Wilry, Berlin , (2000). 

[10] M. O. Scully, M. S. Zubairy, Quantum optics, Cambridge, (1997). 



15 



